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summary 


Nonlinear  flutter  of  a  two-dimensional  airfoil  undergoing  plunging 
and  pitching  motions  is  studied  using  a  time  marching  finite  difference 
scheme.  The  structural  nonlinearity  considered  is  of  the  type  due  to  a  spring 
with  preload  and  freeplay.  Flutter  is  determined  from  solutions  of  the 
structural  dynamic  equations  of  motion  when  either  divergent  or  limited 
amplitude  oscillations  are  encountered.  Case  studies  using  various  airfoil 
parameters  and  values  of  preload  and  freeplay  are  carried  out.  The  effect  of 
initial  condition,  which  is  important  in  nonlinear  problems,  is  investigated 
by  varying  the  displacement  from  equilibrium  of  the  pitch  angle  at  the 
beginning  of  the  airfoil  motion.  For  nonzero  values  of  the  preload,  three 
types  of  oscillatory  motion  are  possible,  namely:  damped,  limited  amplitude 
and  divergent.  The  divergent  flutter  boundary  is  practically  identical  to  that 
for  the  linear  flutter  case.  The  location  of  the  limit-cycle  flutter  boundary 
varies  for  different  airfoil  and  spring  parameters.  For  zero  preload,  damped 
oscillations  are  not  encountered  even  for  air  speeds  down  to  15  percent 
of  the  linear  flutter  speed  which  is  the  lowest  used  in  this  study.  The  limited 
amplitudes  of  the  pitch  and  plunge  motions  are  found  to  be  independent 
of  initial  angular  displacement.  The  characteristics  of  the  oscillations  and 
the  development  of  higher  harmonics  in  the  various  regions  defined  by  the 
flutter  boundary  curves  are  investigated. 


L’aeroelasticite  non-lineaire  d’un  profil  d’aile  bidimensionnel 
subissant  des  mouvements  de  tangage  est  etudiee  au  moyen  d’un  schema  en 
differences  finies  en  fonction  du  temps.  La  non-linearite  structurale  etudiee 
est  celle  produite  par  un  ressort  precharge  et  a  jeu  libre.  L’aeroelasticite 
est  determinee  a  partir  des  solutions  aux  equations  dynamiques  structurales 
de  mouvement  lorsqu’il  y  a  oscillations  a  amplitude  limitee  ou  divergente. 
Des  etudes  de  cas  sont  menees  a  partir  de  divers  parametres  de  profil  d’aile 
et  de  diverses  valeurs  de  precharge  et  de  jeu  libre.  On  etudie  l’effet  de  l’etat 
initial,  important  dans  les  problemes  de  non-linearite,  en  faisant  varier  le 
mouvement  par  rapport  au  point  d’equilibre  de  Tangle  de  tangage  des  que  le 
profil  d’aile  commence  a  se  de  placer.  Pour  des  valeurs  de  precharge  diffe- 
rentes  de  zero,  trois  types  d’oscillation  sont  possibles,  soit:  oscillation 
amortie,  oscillation  a  amplitude  limitee  et  oscillation  divergente.  La  limite 
d’aeroelasticite  divergente  est  pratiquement  identique  a  celle  de  l’aero- 
elasticite  lineaire.  L’emplacement  des  oscillations  limites  d’aeroelasticite 
varie  en  fonction  des  parametres  de  profil  d’aile  et  de  ressort.  Pour  une 
precharge  nulle,  aucune  oscillation  amortie  n’est  relevee,  meme  dans  le  cas 
de  vitesses  aerodynamiques  ramenees  a  15  pour-cent  de  la  vitesse  d’aero¬ 
elasticite  lineaire,  laquelle  est  la  plus  faible  utilisee  dans  le  cadre  de  l’etude. 
On  a  trouve  que  les  amplitudes  limitees  des  mouvements  de  tangage  etaient 
independantes  du  deplacement  angulaire  initial.  Les  caracteristiques  des 
oscillations  et  la  production  d’harmoniques  superieures  dans  les  diverses 
zones  delimitees  par  les  courbes  limites  d’aeroelasticite  sont  etudiees. 


CONTENTS 


SUMMARY 

ILLUSTRATIONS 

SYMBOLS 

TABLE 

APPENDICES 

1.0  Introduction 
2.0  Analysis 

2.1  Two-Degree-of-Freedom  Motion  of  a  2-D  Airfoil 

2.2  Finite  Difference  Scheme 

2.3  Starting  Procedure 

2.4  A  Recurrence  Formula 
3.0  Results  and  Discussions 

3.1  Effect  of  Preload  on  Flutter  Boundary 

3.2  Effect  of  Freeplay  on  Flutter  Boundary 

3.3  Effect  of  Zero  Preload  on  Flutter  Boundary 

3.4  Effect  of  Airfoil-Air  Mass  Ratio  p  on  Flutter  Boundary 

3.5  Effect  of  Uncoupled  Plunge  to  Pitch  Natural 
Frequency  Ratio  ui  on  Flutter  Boundary 

3.5.1  Effect  of  Preload  for  Constant  Freeplay 

3.5.2  Effect  of  Freeplay  for  Constant  Preload 

3.5.3  Effect  of  Airfoil-Air  Mass  Ratio  p 

3.6  Period  of  Limited  Amplitude  Oscillations 
4.0  Conclusions 

References 


Page 

(iii) 

(v) 

(xii) 

29 

73 

1 

3 

3 

6 

8 

9 

10 

11 

15 

17 

18 

19 

19 

20 
21 
21 
22 


(v) 


r.*r s : ' , Tf>tyr  *? 


i 


•  OS1".1' 


5.0 


26 


LIST  OF  ILLUSTRATIONS 


Figure 

1 

2 

3 

4 

5 

6a 

6b 

6c 

6d 

6e 

6f 


6g 

6h 


Two-degree-of-freedom  airfoil  motion 

Nonlinearity  in  pitch  degree-of-freedom  with  preload  and 
freeplay 

Flutter  boundary  for  y  =  100,  u  =  0.2,  hf+  =  0.25°, 
h,-  =  0.75°  and  V,0  =  0.25° 


Flutter  boundary  for  y  =  100,  u  =  0.2,  h^+  =  0.5°,- 
hf .  =1.0°  and  VfQ  =0.5° 


Flutter  boundary  for  y  =  100 

,  u 

=  0 

•  2 ,  h^  t 

hf_  =  1.5°  and 

Vf< 

3  =  1.0“ 

Time  variation 

of 

a  for  a(0) 

= 

8°, 

U*/U*  = 

L 

(corresponding 

to 

location  1 

in 

Fig 

•  3) 

Time  variation 

of 

a  for  a(0) 

= 

8°, 

u*/u*  = 

(corresponding 

to 

location  2 

in 

Fig 

•  3) 

Time  variation 

of 

a  for  a(0) 

= 

8°, 

u*/u*  = 

(corresponding 

to 

location  3 

in 

Fig 

•  3) 

Time  variation 

of 

a  for  a(0) 

= 

8°, 

U*/u*  = 

(corresponding 

to 

location  4 

in 

Fig 

•  3) 

Time  variation 

of 

a  for  a(Q) 

= 

8°, 

U*/U*  = 

L 

(corresponding 

to 

location  5 

in 

Fig 

•  3) 

Time  variation 

of 

a  for  a(0) 

= 

8°, 

u*/u*  = 

L* 

(corresponding 

to 

location  6 

in 

Fig 

•  3) 

Time  variation 

of 

a  for  o(0) 

= 

8°, 

u*/u*  = 

(corresponding 

to 

location  7 

in 

Fig 

•  3) 

Time  variation 

of 

a  for  a( 0) 

= 

-0.3 

°,  U*/U* 

(corresponding 

to 

location  8 

in 

Fig 

•  3) 

Page 

31 

32 

33 

34 

35 

36 

36 

37 

37 

38 

38 

39 


(VI) 


39 


Figure 


Page 


6i  Time  variation  of  a  for  a(0)  =  -1.5°,  U*/U*  =0.78 

(corresponding  to  location  9  in  Fig.  3)  40 

6j  Time  variation  of  a  for  a(0)  =  1.5°,  U*/U*  =  0.84 

Li 

(corresponding  to  location  10  in  Fig.  3)  40 

7a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  u  =  0.2,  hf +  =  0.25°,  hf.  =  0.75°  and 

VfG  =  0.25°  41 

7b  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  w  =  0.2,  hf ♦  =  0.25°,  hf.  =  0.75°  and 

VfC  =  0.25°  41 

8a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  w  =  0.2,  hf «  =  0.5°,  hf .  =  1.0°  and  VfQ  =  0.5°  42 

8b  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  w  =  0.2,  hf ♦  =  0.5°,  hf-  =  1.0°  and  VfD  =  0.5°  42 

9a  Variation  of  limit-amplitude  with  speed  ratio  for 

p  =  100,  w  =  0.2,  hf+  =  1.0°,  hf .  =  1.5°  and  VfD  =  1.0°  43 

9b  Variation  of  limit-amplitude  £  with  speed  ratio  for 

A 

y  =  100,  u  =  0.2,  hf ♦  =  1.0°,  hf.  =  1.5°  and  Vf0  =  1.0°  43 

10  Flutter  boundary  for  y  =  100,  w  =  0.2,  h^,  =  0.5°, 

hf -  =  0.75°  and  Vf0  =  0.5°  44 

11  Flutter  boundary  for  y  =  100,  w  =  0.2,  h^+  =  0.5°, 

hf_  =  1.5°  and  Vf0  =  0.5°  45 

12a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  u  =  0.2,  hf ♦  =  0.5°,  hf.  =  0.75°  and  Vf0  =  0.5°  46 


(vii) 


Figure 


Page 


'Vl 


■«a 

.»*  ‘I 


v: 


'as 


m 


$ 


** 


12b  Variation  of  limit-amplitude  with  speed  ra*io  for 

y  =  100,  u  =  0.2,  hf ♦  =  0.5°,  hf_  =  0.75°  and  V  0  =  0.5°  46 

13a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  u  =  0.2,  hf*  =  0.5°,  hf.  =  1.5°  and  V  «,  =  0.5°  47 

13b  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  u  =  0.2,  hf ♦  =  0.5°,  hf_  =  1.5°  and  Vf0  =  0.5°  47 

14  Flutter  boundary  for  y  =  100,  u>  =  0.2,  h^+  =  0°, 

h, .  =1.0°  and  Vf0  =  0°  48 

15a  Variation  of  limit-amplitude  aA  with  speed  ratio  for 

y  =  100,  u  =  0.2,  hf ♦  =  0°,  h  =  1.0°  and  V  «,  =  0°  49 

15b  Variation  of  limit-amplitude  £  with  speed  ratio  for 

A 

y  =  100,  u  =  0.2,  hf ♦  =  0°,  hf.  =  1.0°  and  VfD  =  0°  49 

16a  Time  variation  of  a  for  a(0)  =  4°,  U*/U*  =  0.16, 

hf*  =  0°,  hf _  =  1.0°  and  vffl  =  0°  50 

16b  Time  variation  of  a  for  a(0)  =  4°,  U*/U£  =  0.636, 

hf.  =  0°,  hf.  =  1.0°  andvf0  =  0°  50 

17  Flutter  boundary  for  y  =  50,  u  =  0.2,  h^»  =  0.5°, 

hf.  =  1.0°  and  V  0  =  0.5°  51 

18  Flutter  boundary  for  y  =  250,  w  =  0.2,  h^»  =  0.5°, 

hf_  =  1.0°  and  V  0  =  0.5°  52 

19a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  50,  u  =  0.2,  hf ♦  =  0.5°,  hf.  =  1.0°  and  V  0  =  0.5°  53 

19b  Variation  of  limit-amplitude  £  with  speed  ratio  for 

A 

y  =  50,  u  =  0.2,  hf ♦  =  0.5°,  h  =  1.0°  and  Vf0  =  0.5°  53 


(viii) 


20a  Variation  of  limit-amplitude  with  speed  ratio  for 
y  =  250,  u  =  0.2,  hf*  =  0.5°,  hf.  =  1.0°  and  Vf0  =  0 

20b  Variation  of  limit-amplitude  with  speed  ratio  for 
y  =  250,  w  =  0.2,  hf*  =  0.5°,  hf.  =  1.0°  and  Vf0  =  0 

21  Flutter  boundary  for  y  =  100,  w  =  0.8,  h^«.  =  0.25°, 
hf .  =  0.75°  and  V<e  =  0.25° 

22  Flutter  boundary  for  y  =  100,  w  =  0.8,  h^+  =  0.5°, 
hf_  =  1.0°  and  Vf0  =  0.5° 

23  Flutter  boundary  for  y  =  100,  w  =  0.8,  h^*  =  1.0°, 
hf .  =  1.5°  and  V  0  =  1-0° 

24a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  u  =  0.8,  hf ♦  =  0.25°,  hf_  =  0.75°  and 

VfC  =  0.25° 

24b  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  u  =  0.8,  hf ♦  =  0.25°,  hf.  =  0.75°  and 

VfQ  =  0.25° 

25a  Variation  of  limit-amplitude  with  speed  ratio  for 
y  =  100,  u  =  0.8,  hf ♦  =  0.5°,  hf.  =  1.0°  and  V  »  =  0 

25b  Variation  of  limit-amplitude  with  speed  ratio  for 
y  =  100,  w  =  0.8,  hf ♦  =  0.5°,  hf.  =  1.0°  and  VfC  =  0 

26a  Variation  of  limit-amplitude  with  speed  ratio  for 
y  =  100,  u  =  0.8,  hf *  =  1.0°,  hf.  =  1.5°  and  VfC  =  1 

26b  Variation  of  limit-amplitude  with  speed  ratio  for 


Figure  Page 

27  Flutter  boundary  for  y  =  100,  u  =  0.8,  =  0.5°, 

hf-  =  0.75°  and  VfC  =0.5°  61 

28  Flutter  boundary  for  y  =  100,  u  =  0.8,  =  0.5°, 

hf.  =  1.5°  and  Vf0  =  0.5°  62 

29a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  u  =  0.8,  hf ♦  =  0.5°,  hf.  =  0.75°  and  VfG  =  0.5°  63 

29b  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  w  =  0.8,  hf ♦  =  0.5°,  hf.  =  0.75°  and  Vf0  =  0.5°  63 

30a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  100,  w  =  0.8,  hf ♦  =  0.5°,  hf.  =  1.5°  and  VfC  =  0.5°  64 

30b  Variation  of  limit-amplitude  £  with  speed  ratio  for 

y  =  100,  u  =  0.8,  hf ♦  =  0.5°,  hf.  =  1.5°  and  Vf0  =  0.5°  64 

31  Flutter  boundary  for  y  =  50,  w  =  0.8,  h^*  =  0.5°, 

hf.  =  1.0°  and  VfD  =  0.5°  65 

32  Flutter  boundary  for  y  =  250,  u  =  0.8,  h^*  =  0.5°, 

hf_  =  1.0°  and  VfC  =  0.5°  66 

33a  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  50,  u  =  0.8,  hf *  =  0.5°,  hf.  =  1.0°  .ad  VfQ  =  0.5°  67 

33b  Variation  of  limit-amplitude  with  speed  ratio  for 

y  =  50,  u  =  0.8,  hf *  =  0.5°,  hf.  =  1.0°  and  VfQ  =  0.5°  67 

34a  Variation  of  limit-amplitude  a with  speed  ratio  for 

y  =  250,  u  =  0.8,  hf*  =  0.5°,  hf-  =  1.0°  and  Vf0  =  0.5°  68 

34b  Variation  of  limit-amplitude  £  with  speed  ratio  for 

y  =  250,  u  =  0.8,  hf  ♦  =  0.5°,  hf.  =  1.0°  and  VfD  =  0.5°  68 


(x) 


Figure 


Page 


35  Period  of  limited  amplitude  oscillations  versus  speed 

ratio  for  p  =  100,  w  =  0.2  69 

36  Period  of  limited  amplitude  oscillations  versus  speed 

ratio  for  p  =  100,  u  =  0.8  70 

37  Period  of  limited  amplitude  oscillations  versus  speed 

ratio  for  p  =  100,  u  =  0.2  71 

Appendices 

Appendix  Page 

A  73 

B  77 

Table 

Table 

1  List  of  Airfoil  and  Nonlinear  Spring  Parameters  used  in 

Case  Studies  29 


i 


» 

» 


1 

n 

N 

X 


IPu 


Symbol 


Definition 


constant  in  Equation  (11) 

nondimens ional  distance  measured  from  airfoil 
mid-chord  to  elastic  axis 

semi-chord  of  airfoil;  also  constant  in  Equation  (11) 
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1.0  Introduction 


The  assumption  of  structural  linearity  is  frequently  made  to 
enable  the  use  of  analytical  methods  in  determining  the  divergence  and 
flutter  characteristics  of  aerodynamic  surfaces.  Linear  theory  will 
predict  the  magnitude  of  dynamic  pressure  or  flight  velocity  above 
which  the  system  under  consideration  becomes  unstable  and  the  motion 
grows  exponentially  in  time.  However,  aircraft  structures  often 
exhibit  nonlinearities  which  affect  not  only  the  flutter  speed,  but 
also  the  characteristics  of  flutter  motion.  An  understanding  of  the 
nonlinear  behaviour  of  the  system  is  crucial  to  the  efficient  and  safe 
design  of  aircraft  wings  and  control  surfaces. 

Reference  1  gives  an  excellent  discussion  of  the  various 
types  of  nonlinearities  and  some  of  the  methods  used  to  treat  them.  In 
general,  structural  nonlinearities  can  be  categorized  as  either 
distributed  or  concentrated.  Only  the  latter  is  considered  in  this 
study  because  they  are  more  important.  An  example  of  this  type  of 
nonlinearity,  which  is  investigated  in  this  report,  is  the  freeplay 
representing  a  loose  hinge  or  linkage  of  a  control  system.  The 
response  characteristics  are  usually  functions  of  the  amplitude  of 
oscillation.  The  system  stiffness  behaviour  changes  with  amplitude  of 
motion  and  at  some  particular  flight  speed  the  motion  can  be 
self-excited  and  attain  a  limited  amplitude.  The  appearance  of  the 
phenomenon  of  limit-cycle  flutter  is  important  from  the  design 
viewpoint.  These  oscillations  may  occur  within  the  divergence  and 
flutter  flight  envelope,  but  the  amplitude,  frequency  and  duration  of 
these  limit-cycle  oscillations  may  have  an  important  impact  on  the 
structural  integrity  of  the  aerodynamic  surfaces. 


Woolston  et  al  (Ref.  2)  studied  the  effects  of  structural 


nonlinearities  on  the  flutter  of  a  two-dimensional  airfoil  using  an 
analog  computer.  There  are  serious  drawbacks  in  the  use  of  an  analog 
computer  to  analyse  nonlinear  flutter  and  accuracy  is  often  not  as 
great  as  one  would  desire.  An  alternate  analytical  approach  was 
suggested  by  Shen  (Ref.  3)  using  the  well  known  Kryloff  and  Bogoliuboff 
(Ref.  4)  method  in  nonlinear  vibration  theory.  The  original  limitation 
of  weak  nonlinearities  was  removed  by  adopting  a  modification  given  by 
Popov  (Ref.  5).  Only  rigid  surfaces  with  single  nonlinearity  were 
treated.  Extensions  to  multiple  nonlinearities  (Ref.  6,  7)  and 
inclusion  of  higher  harmonic  terms  (Ref.  8)  have  been  reported  in 
recent  studies. 

Another  method  of  investigating  nonlinear  flutter  is  by 
numerical  time  integration  of  the  structural  equations  of  motion.  This 
has  the  advantage  that  no  assumptions  on  the  type  of  motion  have  to  be 
made  beforehand.  Also,  the  effect  of  initial  conditions,  which  is 
important  in  nonlinear  problems,  can  readily  be  investigated. 

Flutter  analyses  of  two-dimensional  airfoils  with  aerodynamic 
nonlinearities  and  structural  nonlinearities  of  the  type  due  to  a  cubic 
spring  have  been  studied  by  Lee  and  LeBlanc  (Refs.  9  and  10).  In  that 
study,  Houbolt's  (Ref.  11)  finite  difference  scheme  was  found  to  be 
efficient  and  of  sufficient  accuracy  in  determining  flutter 
boundaries.  Following  the  approach  used  in  Ref.  10,  the  aeroelastic 
equations  for  a  two-dimensional  airfoil  performing  plunging  and 
pitching  motions  are  written  in  this  report  in  finife  difference  form 
using  incompressible  aerodynamics  in  the  time  domain  given  by  Wagner's 
function.  The  effect  of  initial  conditions  on  nonlinear  flutter  is 


studied  by  varying  the  displacement  from  equilibrium  of  the  pitch  angle 
at  the  beginning  of  the  airfoil  motion.  Divergent  and  limit-cycle 
flutter  are  investigated  using  different  airfoil  parameters  and  values 
of  spring  preload  and  freeplay.  The  characteristics  of  the  oscilla¬ 
tions  and  the  development  of  higher  harmonics  in  the  various  regions 
defined  by  the  flutter  boundary  curves  are  analysed.  The  results  are 
useful  not  only  in  assessing  the  validity  and  accuracies  of  the  method 
of  harmonic  balance  (Refs.  3-8)  commonly  used  in  the  study  of  nonlinear 
flutter,  but  also  in  understanding  the  impact  of  limit-cycle 
oscillations  on  the  structural  integrity  of  aerodynamic  surfaces. 

2 . 0  Analysis 

2 . 1  Two-Degree-of-Freedom  Motion  of  a  2-D  Airfoil 

Figure  1  shows  the  notations  used  in  the  analysis  of  a 
two-degree-of-f reedom  motion  of  an  airfoil  oscillating  in  pitch  and  in 
plunge.  The  plunging  deflection  is  denoted  by  h,  positive  in  the 
downward  direction,  a  is  the  pitch  angle  about  the  elastic  axis, 
positive  with  the  nose  up.  The  elastic  axis  is  located  at  a  distance 
a^b  from  the  midchord,  while  the  mass  centre  is  located  at  a  distance 
x^b  from  the  elastic  axis.  Both  distances  are  positive  when  measured 
towards  the  trailing  edge  of  the  airfoil.  The  aeroelastic  equations  of 
motion  for  linear  springs  have  been  derived  by  Fung  (Ref.  12).  For 
nonlinear  restoring  moment  from  a  spring  with  preload  and  freeplay  in 
the  torsional  degree  of  freedom  (Figure  2) ,  they  can  be  written  as 


fol lows : 
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mh  +  Sa  ♦  C.  h  + 
n 

Sh+Ia+Ca+  F(a)  =  R  (2) 

a  a  v  ' 

where  the  dot  represents  differentiation  with  respect  to  time.  S  is 

the  airfoil  static  moment  about  the  elastic  axis,  F(a)  is  a  nonlinear 

function  representing  the  restoring  moment,  P  and  R  are  the  externally 

applied  force  and  moment  respectively.  The  symbols  m,  C^,  K^,  I  and 

C&  are  the  airfoil  mass,  linear  damping  coefficient,  linear  spring 

constant,  wing  mass  moment  of  inertia  about  the  elastic  axis,  and 

torsional  damping  constant  respectively. 

Define  £  =  h/b,  X(i  »  S/bm,  ^  =  (Kh/m)*,  u>a  =  (K^)*, 

r  =  (I  /mb*)*,  =  C,/2(mK.  )*,  C  =  C  /2(I  K  )*.  Here  K  is  the 

aa  t,  h  n  a  a  a  a  a 

torsional  spring  constant  for  the  linear  part  of  the  moment- 
displacement  curve  shown  in  Figure  2.  Equations  (1)  and  (2)  can  be 
written  in  nondimensional  form  as  follows: 


V 


=  p 


(1) 


5"(T)  +  xaa"(T)  +  2^  gs  r(t)  +  [  —  }'  5(T)  =  P(T)  (3) 

A  r(T)  +  a"(T)  +  2C  ^lill  +  J-  F[a(t)  ]  =  r(t)  (4) 

r  U* 

a 

where  p(t)  and  r(t)  are  the  nondimensional  force  and  moment 
respectively.  The  prime  denotes  differentiation  with  respect  to  the 
non-dimensional  time  t  defined  as 
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T  = 


Ut 


(5) 


In  Eqs.  (3)  and  (4),  u  is  given  by 


U*  is  defined  as 


-  ^ 
U)  =  — 1 

u 

a 


(6) 


U*  = 


U 

bu 


(7) 


From  Figure  2.,  F(a)  can  be  written  as 


F(«)  = 


a  +  V,0  -  h,»  a  <  h,. 


Vfo 


+  Vf,  -  hf. 


h^.<  a 


(8) 


For  incompressible  flow,  Fung  (Ref.  12)  gives  the  following 
expressions  for  p(t)  and  r(x): 


P(*)  =  ~  ^[?"(t)-aha"(t)+a'  (t)] 


-  n[fa(0)+C' (O)f(i-a.  )a'(0)}*(T) 


+  /  *(T-o)[a’(o)+C(o)+(i-a  )a"(a)]do] 
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r(T)  =  -J-y(Wh){[a(0)+r  (0)  +  (i-a  )a'(0)]f(T) 

M  a  11 

j 

+/  $(x-a) [a' (o)+£"(a)+($-ah)a"(a) ] da } 
o 

+  ^{vr(T)~aha’,{T)1~(i~ah)a,(T)  ~  ia”(T))  (10) 

^  a 


where  p  is  the  airfoil-air  mass  ratio,  and  the  Wagner's  function  $(x) 
is  given  by 


$(t)  =  £  -  ae  -  ce  (11) 


and  for  incompressible  flow,  the  constants  £  =  1,  a  =  0.165,  b  = 
0.0455,  c  =  0.335  and  d  =  0.3  are  given  in  Ref.  12. 


2 . 2  Finite  Difference  Scheme 

For  a  system  with  few  structural  components  or  the  number  of 
vibration  modes  is  small,  higher  order  methods  such  as  the  eighth  order 
scheme  reported  in  Ref.  13  does  not  offer  any  distinct  advantage  over 
Houbolt's  (Ref.  11)  scheme  which  is  simpler  and  less  cumbersome  to  use. 
For  this  reason,  Houbolt's  method  is  used  in  the  present  analysis.  The 
derivatives  at  time  x+Ax  are  replaced  with  backward  difference  formulas 
using  values  at  three  previous  points.  For  example, 


<x*  (t+At )  =  ^  [  1  la( t+At )-18a( t )+9a(t-At )-2a(t-2At )  ]  (13) 


Similar  expressions  can  be  written  for  £"(t+At)  and  £'(t+At).  In 
difference  form,  Eqs .  (3)  and  (4)  can  be  expressed  as: 


Pna(t+At)  +  Pl2  £(t+At)  =  X1 


(14) 


[P21  +  U*  2  F  ] a( t  +Ax )  +  P22£(i+At)  =  j^-U*  2Fa  (15) 


In  Eq.  (15),  F[a(t)]  has  been  replaced  with  F^(a;a)a(T+AT ) 


+  F  (a) ,  where 
a 


Fp(a;a)  = 


a  <  hf  ,  ,  hf .  <  a 


otherwise 


(16) 


and 
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Renaming  the  coefficients  in  Eqs .  (14)  and  (15)  yields  the 

following : 


{ a(t+At ) 

X1 

[P]  { 

.  =  . 

|  £(t+At ) 

ix2) 

(18) 


which  has  for  solution 


a(t+At) 

£(t+At) 


(19) 


[P]  can  take  on  one  of  the  two  forms  whereas  (X)  can  take  on  one  of 
three  forms  depending  on  the  estimated  value  of  a(t+Ai). 


2 . 3  Starting  Procedure 


Houbolt's  scheme  requires  values  of  a  and  £  at  times  t-2At , 
t-At  and  t  in  order  to  determine  their  values  at  t+At .  At  time  t  =  0, 
a  special  starting  procedure  is  required.  A  Taylor  series  is  used  to 
obtain  the  following: 


a(-At)  =  a(0)  -  Atcx'(O)  +  j-  a"(0)  +  0(At3)  (20) 

A  2  3 

a( At )  =  a(0)  +  Aia'(O)  +  a"(0)  +  0(At  )  (21) 


with  similar  expressions  for  £(-At)  and  £(At).  The  expressions  for  X 
and  in  Appendix  A  have  terms  depending  on  the  first  and  second 
derivatives  of  a  and  £,  and  these  derivatives  are  given  as  follows: 


1 


a’(-At)  =  a' (0)  -  At  a" ( 0 )  +  a"(0)  +  0 ( At  3 ) 


(22) 
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a' (Ax)  =  a'(0)  +  Axa"(0)  +  j-  a"’(0)  +  0(Ax3)  (23) 


and 


a"(-AT)  =  a" (0)  -  At  a'”(0)  +  j-  aIV(0)  +  0(At3)  (24) 

2 

a"(Ax)  =  a"(0)  +  Ax  a'"(0)  +  cxIV(0)  +  0(Ax3)  (25) 


with  expressions  of  the  £  derivatives  written  in  the  same  form  as  those 
given  in  Eqs .  (22)  to  (25).  The  initial  conditions  a(0) ,  a'(0),  £(0) 
and  £'(0)  are  known,  and  the  higher  derivatives  can  be  obtained  from 
Eqs.  (3)  and  (4)  and  are  given  by: 


IQ] 


c(n)(0)l 

Y  (n)' 
1 

i(n)(0) j 

Y  (n) 

1 2 

which  has  for  solution 


r.(n,(o,' 

•  =  [Qf1  • 

IV-j 

k<n,j 

(26) 


(27) 


Here  n  =  2,  3  and  4.  Q^,  Q^,  Q^,  and  Y^^  are  given  in 

Appendix  B.  For  the  next  step,  Houbolt's  scheme  can  be  used  since  a 
and  £  at  t  =  -Ax,  0  and  Ax  are  known. 


2 . 4  A  Recurrence  Formula 

In  Eq.  (18),  the  terms  and  X2  contain  integrals  which  have 
to  be  evaluated  at  each  time  step.  To  reduce  the  amount  of 


10 


computations,  a  recurrence  formula  is  used.  Using  Simpson's  rule,  the 
integral  I^(x+Ax)  in  Eq.  (A6)  of  Appendix  A  can  be  written  as 

I1(t+At)  =  e-bATI1(i)  +  (9X(t+At)  +  19X(x)e“bAT 

-  5X(T-At)e"2bAT  +  X(T-2At)e'3bAl}  (28) 

where  X  is  defined  in  Eq.  (A8) . 

The  starting  value  for  I^(Ax)  is  given  by 

I x ( At )  =  ~  { 2 X ( At )  +  5X(0)e"bAT  -  X(-Ax)e~2bAt }  (29) 

Expressions  for  ^(t+Ax)  and  I2(Ax)  are  obtained  by  substituting  b  for 
d  in  the  above  equations. 

3 . 0  Results  and  Discussions 

In  this  report,  only  nonlinearities  in  the  pitch  degree  of 

freedom  are  considered.  The  elastic  axis  of  the  airfoil  is  placed  at 

the  j  chord  location  (that  is,  a,  =  -0.5),  r  and  x  are  kept  constant 

h  a  a  r 

at  0.5  and  0.25  respectively.  The  two  properties  of  the  airfoil  being 
varied  are  p  and  u>.  The  effects  of  preload  and  freeplay  are 
investigated  by  varying  the  values  of  V^0,  h^*  and  h^..  Table  1  shows 
the  combinations  of  parameters  used  in  this  study. 

A  comparison  of  cases  1,  2  and  3  shows  the  effect  of  preload 
for  constant  freeplay,  while  cases  2,  4  and  5  show  the  effect  of 
freeplay  for  constant  preload.  Cases  6  and  7  are  with  zero  preload  and 
the  effect  of  varying  p  can  be  obtained  from  cases  2,  8  and  9.  The 
results  for  increasing  w  to  a  value  of  0.8  are  given  in  cases  10  to  16 


where  the  values  of  preload  and  freeplay  are  kept  the  same  as  those  for 
cases  1  to  5,  8  and  9  in  which  u  =  0.2.  These  computations  are  carried 
out  primarily  to  study  the  effect  of  u  on  the  binary  flutter  of  the 
airfoil . 

3 . 1  Effect  of  Preload  on  Flutter  Boundary 

Figures  3  to  5  show  the  flutter  boundaries  for  cases  1  to  3. 
The  values  of  y  and  w  are  kept  constant  at  100  and  0.2  respectively. 

The  freeplay  is  0.5°  with  the  preload  V^0  set  at  0.25°,  0.5°  and  1° 
respectively.  Throughout  this  study,  it  is  assumed  that  h^*  =  V^0  so 
that  the  restoring  moment  is  zero  when  the  displacement  a  =  0. 

To  determine  the  flutter  boundary,  Eq.  (19)  is  solved  for 
given  initial  conditions.  In  this  report,  only  a(0)  is  varied  and 
£(0)  =  S'(0)  =  a'(0)  =  0.  The  procedure  is  to  find  the  linear  flutter 
speed  U*  first,  and  this  is  equivalent  to  solving  the  problem  for 
h^_  =  h^*.  In  the  nonlinear  case,  once  a(0)  is  specified,  a  value  of 
U*  usually  greater  than  U*  is  selected  and  a  and  £  are  obtained  by  the 
time  marching  finite  difference  scheme.  For  the  type  of  structural 
nonlinearity  considered,  the  solution  is  divergent  for  U*  >  U*  and  the 
nonlinear  divergent  flutter  speed  is,  within  the  numerical  accuracy, 
the  same  as  U*  for  these  three  cases  considered.  It  is  observed  that 
for  a(0)  >  V  o  decreasing  U*  below  U*  results  in  limit-cycle  flutter. 
The  oscillation  is  sel f-excitated  and  maintains  a  constant  amplitude 
which  is  self-limited.  The  boundaries  between  divergent  and 
limit-cycle  flutter  in  Figures  3  to  5  are  practically  vertical  lines  at 
U*/U*  =  1  for  values  of  a(0)  ranging  from  -10°  to  20°. 

As  U*  decreases,  a  value  will  be  reached  where  any  further 


decrease  will  result  in  damped  oscillations  of  the  airfoil.  Boundaries 
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can  be  identified  on  the  a(0)  versus  U*/U£  plots  separating  the  regions 
of  limit-cycle  flutter  with  the  stable  regions  where  the  oscillations 
decay . 

Figures  3  to  5  also  show  that  for  a  freeplay  of  0.5°  the 
region  of  limit-cycle  flutter  decreases  with  increasing  values  of  the 
preload.  Another  observation  is  the  non-symmetry  of  the  boundaries 
between  decaying  oscillation  and  limit-cycle  flutter  with  initial 
displacement  a(0) .  For  values  of  initial  displacement  less  than  the 
preload,  the  system  moves  on  the  linear  part  of  the  moment- 
displacement  curve  shown  in  Figure  2,  and  at  velocities  below  the 
linear  flutter  speed,  the  system  is  stable.  Increasing  a(0)  will  have 
a  destablizing  effect,  but  this  is  only  restricted  to  the  limited 
amplitude  oscillation  regions. 

Figure  3  shows  pockets  in  the  damped  oscillation  region  where 
the  airfoil  oscillates  with  constant  amplitude.  This  limit-cycle 
flutter  behaviour  is  only  observed  for  the  smallest  preload  value  of 
0.25°  and  is  not  detected  in  the  other  two  cases.  The  regions  of 
limited  amplitude  oscillations  are  determined  using  a  binary  search 
complemented  with  linear  grid  scans.  This  by  no  means  assures  that  all 
such  regions  especially  the  small  ones  have  been  identified,  but  those 
that  are  found  can  be  considered  to  be  quite  accurately  defined. 

Figures  6a  to  6j  show  the  time  behaviour  of  a  in  the  various 
regions  corresponding  to  the  locations  marked  1  to  10  in  Figure  3. 
Figure  6a  is  for  initial  displacement  a(0)  =  8°  and  U*  =  0.7U*.  The 
airfoil  damps  out  fairly  rapidly  and  the  oscillations  are  mainly  sinu¬ 
soidal  and  no  noticeable  harmonics  are  detected.  This  is  a  typical  a 
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decay  curve  sufficiently  far  from  any  limit-cycle  futter  region.  The 
same  behaviour  is  also  observed  for  the  £  oscillations. 

Figure  6b  corresponds  to  location  2  in  Figure  3  where 
a(0)  =  8°  and  U*/U*  =  0.78.  This  type  of  oscillatory  behaviour  of  a  is 
typical  in  those  pockets  where  limit-cycle  flutter  occurs.  A  strong 
second  harmonic  is  always  present. 

At  initial  displacment  a(0)  =  8°  and  U*/U*  =  0.81, 
corresponding  to  location  3  in  Figure  3,  the  a  decay  curve  is  shown  in 
Figure  6c.  Higher  harmonics  can  be  detected  from  the  oscillatory 
motion  of  the  airfoil,  and  this  is  often  the  case  in  the  proximity  of 
the  limit-cycle  flutter  regions.  Figure  6d  is  also  in  the  damped 
oscillation  region  but  it  is  very  close  to  the  limit-cycle  flutter 
boundary.  This  corresponds  to  location  4  in  Figure  3  with  initial 
displacement  a(0)  =  8°  and  U*/U£  =  0.82.  The  airfoil  oscillates  in 
pitch  (also  in  plunge)  at  approximately  constant  amplitude  for  a  number 
of  cycles  (nearly  twenty  in  this  case)  and  suddenly  damps  out.  This  is 
often  encountered  very  close  to  the  limit-cycle  flutter  boundary. 

Figure  6e  shows  the  time  variation  of  a  in  the  limit-cycle 

flutter  region  for  a(0)  =  8°  and  U*/U*  =  0.S3  corresponding  to  location 

5  in  Figure  3.  Again,  being  close  to  the  flutter  boundary,  the 

appearance  of  a  second  harmonic  can  be  detected.  Sufficiently  far  away 

from  the  limit-cycle  flutter  boundary,  the  curves  shown  in  Figures  6f 

and  6g  for  limited  amplitude  oscillations  corresponding  to  locations  6 

and  7  in  Figure  3  indicate  the  motion  of  the  airfoil  to  be  practically 

sinusoidal.  This  type  of  motion  is  typical  even  if  U*/U*  is  close  to 

JLi 

the  divergent  flutter  boundary. 


Figure  6h  shows  the  a  oscillations  at  location  8  on  Figure  3 
for  a(0)  =  -0.3°  and  U*/U*  =  0.99785.  For  this  value  of  initial 
displacement,  the  airfoil  motion  is  on  the  linear  part  of  the  moment- 
displacement  curve  shown  in  Figure  2.  This  type  of  decay  curve  is 
typical  for  linear  flutter  (Ref.  10). 

The  behaviour  of  the  airfoil  motion  inside  the  limit-cycle 
flutter  region  at  location  9  in  Figure  3  is  given  in  Figure  6i.  The 
values  of  a(0)  and  U*/U*  are  -1.5°  and  0.78  respectively.  The 
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oscillations  behave  in  a  similar  manner  as  those  given  in  Figure  6b. 
This  is  characterized  by  the  presence  of  a  strong  second  harmonic 
component.  It  is  noted  that  in  all  pockets  of  limnit-cycle  flutter 
within  the  damped  oscillation  region,  the  oscillations  of  the  airfoil 
are  similar  to  those  given  in  Figure  6b  and  6i . 

It  is  observed  in  Figure  3  that  there  are  bulges  on  the 
limit-cycle  flutter  boundary.  For  example,  for  positive  values  of  a(0) 
and  0.825  <  U*/U*  <  0.86,  the  oscillations  can  either  be  of  limited 
amplitude  or  damped  depending  on  the  value  of  a(0).  Figure  6j  shows 
the  time  variation  of  a  corresponding  to  location  10  in  Figure  3  where 
a(0)  =  1.5  and  U*/U*  =  0.84. 

In  Figures  7-9,  the  limited  amplitudes  and  are  plotted 
against  the  speed  ratio  U*/b£  for  cases  1-3.  It  is  found  that  they  are 
independent  of  the  initial  displacement  a(0) .  In  Figure  7,  the  curves 
have  breaks  at  U*/U*  approximately  equal  to  0.825.  To  the  left  of  this 
value  of  U*/U*,  the  amplitudes  are  obtained  from  the  limit-  cycle 
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flutter  pockets  inside  the  damped  oscillation  region.  The  curves  do 
not  join  smoothly  with  those  for  U*/U?  greater  than  0.825  and  this 

Lt 


5 


is  due  to  the  different  oscillatory  motion  of  the  airfoil  as  shown  in 
Figures  6b  and  6f.  A  comparison  of  the  amplitudes  for  U*/U*  >  0.87  in 
Figure  7  with  those  in  Figures  8  shows  that  they  are  identical. 
Similarly,  the  results  in  Figure  9  coincide  with  those  in  Figure  7  for 
U*/U*  >  0.915.  This  indicates  that  the  variations  of  the  limited 
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amplitudes  a  and  £  with  speed  ratio  for  various  values  of  preload, 
but  with  the  freeplay  fixed,  lie  on  the  same  curves.  The  values  of 
U*/U*  where  the  curves  begin  depend  on  the  values  of  the  preload. 

Throughout  this  study,  the  value  of  the  time  step  At  is  taken 
to  be  1/128  of  the  shorter  period  of  the  two  coupled  modes  of 
oscillation  of  the  airfoil  in  the  absence  of  aerodynamic  forces.  The 
periods  of  the  coupled  modes  are  obtained  from  Equation  (3)  and  (4)  for 
F[a(t)]  =  a(t),  and  p(t)  =  r(t)  =  0.  In  numerical  time  integration 
schemes,  Bathe  and  Wilson  (Ref.  14)  pointed  out  that  the  amplitude 
decays  due  to  numerical  errors  and  is  dependent  on  At .  The  value  of  At 
used  in  this  study  is  found  to  be  sufficiently  small  to  give  good 
accuracies  in  determining  the  flutter  boundaries  while  ensuring  the 
computation  time  is  not  excessive. 

3 . 2  Effect  of  Freeplay  on  Flutter  Boundary 

To  investigate  the  effect  of  freeplay  on  the  flutter 
boundary,  results  for  cases  2,  4  and  5  of  Table  1  are  shown  in  Figures 
4,  10  and  11.  The  preload  is  constant  at  0.5°,  and  the  values  of  the 
freeplay  are  0.25°,  0.5°  and  1°  respectively.  The  boundaries  between 
divergent  and  limit-cycle  flutter  are  all  located  at  U*/U*  -  1. 
Decreasing  the  freeplay  will  decrease  the  limited  amplitude 
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oscillation  region  and  move  the  limit-cycle  flutter  boundary  closer  to 

the  divergent  flutter  boundary.  In  the  limit,  as  the  freeplay  tends  to 

zero,  the  two  flutter  boundaries  coincide  at  U*/U*  =  1  which  is  to  be 

Li 

expected  since  the  system  is  acted  on  by  linear  spring  forces.  Because 
of  the  preload,  the  flutter  boundaries  are  not  symmetrical  about 
a (0 )  =  0. 

The  limited  amplitudes  a ^  and  motions  are  given  in  Figures 

8,  12  and  13.  They  are  independent  of  the  initial  displacement  a(0) , 

but  for  given  speed  ratio,  higher  amplitudes  are  obtained  when  the 

values  of  the  freeplay  are  increased.  The  curves  in  Figure  13  have 

breaks  at  U*/U*  approximately  equal  to  0.825.  The  amplitudes  for  U*/U* 
Li  L* 

<  0.825  are  determined  from  the  pockets  of  limit-  cycle  flutter  in  the 
damped  oscillation  region.  The  presence  of  a  strong  second  harmonic  is 
always  observed  in  the  oscillatory  motion  of  the  airfoil  inside  these 
pockets  (see  Figs.  6b  and  6i) .  This  is  quite  unlike  the  motions  for 
(j*/U*  >  0.825  where  the  oscillations  are  mainly  sinusoidal. 

It  is  interesting  to  note  that  when  Figure  3  is  compared  with 
Figure  11,  the  flutter  boundaries  are  identical  if  the  vertical  scale 
for  a(0)  in  Figure  3  is  multiplied  by  a  factor  of  2  which  is  the  ratio 
of  the  two  values  of  preload.  Similarly,  the  limited  amplitudes  and 

£A  in  Figure  7  when  multiplied  by  2  match  those  given  in  Figure  13. 

From  the  limited  results  of  the  two  cases  considered  here,  it  appears 
that  for  a  particular  combination  of  preload  and  freeplay,  similar 
results  are  obtained  for  other  values  of  preload  and  freeplay  if  the 
ratio  of  preload  to  freeplay  is  kept  the  same.  The  flutter  boundary 
curves  can  be  made  identical  if  the  vertical  scales  are 
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multiplied  by  factors  equal  to  the  ratio  of  the  preloads.  Also,  the 
limited  amplitudes  of  the  plunge  and  pitch  oscillations  differ  from 
those  with  other  values  of  preload  and  free  play  by  the  same  factors. 

3 . 3  Effect  of  Zero  Preload  on  Flutter  Boundary 

For  zero  preload,  cases  6  and  7  of  Table  1  are  investigated 
for  a  freeplay  of  1°.  The  difference  bewtween  these  two  cases  is  in 
the  location  of  the  equilibrium  position.  Case  6  gives  the  equilibrium 
of  the  airfoil  at  a  point  on  the  linear  portion  of  the  moment- 
displacement  curve  (Figure  2)  just  ahead  of  the  freeplay,  while  in  case 
7,  the  equilibrium  position  is  centered  in  the  freeplay. 

The  flutter  boundary  for  case  6  is  given  in  Figure  14.  The 
boundary  originates  at  U*/U*  =  1  and  a(0)  =  0  and  rapidly  becomes  a 
straight  line  at  U*/U?  =  0.995.  Unlike  the  previous  cases  with 

la 

positive  preload,  the  figure  shows  only  two  regions,  namely;  divergent 

and  limit-cycle  flutter.  In  this  study,  no  numerical  computations  have 

been  carried  out  for  U*/U*  <  0.15  since  the  number  of  cycles  needed  to 

L 

determine  the  limited  amplitude  becomes  larger  and  larger  as  the  speed 
ratio  decreases,  and  no  drastic  change  in  the  behaviour  of  the  airfoil 
motion  is  anticipated. 

The  limited  amplitudes  a  and  £  are  plotted  in  Figure  15. 

A  A 

The  curves  are  discontinuous  at  U*/U*  =  0.69  and  Figure  16  show  typical 

la 

pitch  oscillations  of  the  airfoil  near  the  discontinuity  (Figure  16b) 
and  at  a  fair  distance  away  (Figure  16a) .  The  mean  value  of  a  is  not  a 
constant  but  depends  on  U*/U*,  and  the  oscillatory  motion  of  the 
airfoil  is  non-symmetrical  about  the  mean.  For  U*/U*  slightly  greater 
than  0.69,  a  small  second  harmonic  can  be  detected,  but  further 
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increase  in  U*/U*  results  in  oscillations  which  are  practically 
sinusoidal.  For  U*/U*  >  0.69,  the  mean  value  of  a  is  a  constant  0.5° 

Li 

which  is  half  the  value  of  the  freeplay. 

For  zero  preload  and  equilibrium  position  centered  in  the 
freeplay,  the  flutter  boundary  and  limited  amplitude  curves  are 
identical  to  Figures  14  and  15  illustrating  their  dependence  only  on 
the  value  of  the  freeplay.  The  mean  motion  for  U* / U*  <  0.69  is  usually 
slightly  above  a  =  0  although  for  some  values  of  U*/U*  and  a(0) ,  it  can 
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be  slightly  below  a  =  0.  However,  for  U*/U*  >  0.69  the  mean  value  of  a 
is  always  zero.  Similar  to  case  6,  a  large  second  harmonic  is  present 
in  the  oscillatory  motion  of  the  airfoil  for  U*/U*  <  0.69.  The 
magnitude  of  the  second  harmonic  is  less  pronounced  in  the  vicinity  of 
y-.'c/jj*  =  0.69  and  as  U*/U*  increases,  its  value  diminishes.  The 
oscillations  appear  mainly  sinusoidal  with  no  noticeable  harmonics  when 
U*/U*  >0.8. 

3 . 4  Effect  of  Airfoil-Air  Mass  Ratio  p  on  Flutter  Boundary 

The  effect  of  airfoil-air  mass  ratio  on  the  flutter  boundary 
for  p  =  50,  100  and  250  are  shown  in  Figures  17,  4  and  18  respectively. 
The  preload  and  freeplay  are  both  0.5°  in  all  three  cases. 

The  boundaries  for  divergent  flutter  are  all  approximately  at 
U*/u*  =  1  for  all  three  values  for  p,  while  those  for  limit-cycle  flut¬ 
ter  move  towards  U*/U;f  -  1  for  increasing  p.  At  p  =  250,  pockets  of 

1-* 

limited  amplitude  oscillations  appear  in  the  damped  oscillation  region. 
The  oscillatory  motions  of  the  airfoil  in  the  various  regions  on  the 
flutter  boundary  curves  are  similar  to  those  given  in  Figure  6. 
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The  limited  amplitudes  aA  and  are  given  in  Figures  19,  8 


and  20  for  the  three  values  of  p  considered.  Comparison  of  the  curves 


for  a shows  that  the  limited  amplitude  between  the  limit-cycle  and 


divergent  flutter  boundaries  practically  does  not  change  with  p  and  no 


differences  can  be  detected.  On  the  other  hand,  £  is  found  to 
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increase  with  p.  This  behaviour  where  p  has  an  effect  on  the  plunge 
degree  of  freedom  but  with  no  noticeable  effect  on  the  pitch  motion  is 
also  observed  for  a  cubic  nonlinearity  in  the  restoring  moment  reported 
in  Ref.  10. 


3 . 5  Effect  of  Uncoupled  Plunge  to  Pitch  Natural  Frequency  Ratio  ui 
on  Flutter  Boundary 


3.5.1  Effect  of  Preload  for  Constant  Freeplay 
In  cases  10  to  12  given  in  Table  1,  the  value  of  u  is  0.8  and 
the  preload  is  set  at  0.25°,  0.5°  and  1°  for  a  constant  freeplay  of 
0.5°.  The  flutter  boundaries  are  shown  in  Figures  21  to  23,  and 
comparison  with  Figures  3  to  5  for  the  corresponding  values  of  preload 
and  freeplay  at  u)  =  0.2  gives  the  effect  of  u>  on  the  oscillatory  motion 
of  the  airfoil. 

From  Figures  3  and  21,  bringing  the  plunge  and  pitch 
uncoupled  natural  frequencies  closer  to  each  other  seem  to  eliminate 
the  porkets  where  limit-cycle  flutter  occurs.  The  bulges  on  the  limit- 
'  V'  it'  flutter  boundary  grow  with  the  preload  for  w  =  0.8.  Similar  to 
i  ';<•  observations  from  cases  1  to  3,  it  is  found  that  the  limited 


impiitudes  a  and  £  given  in  Figures  24-26  are  independent  on  initial 
A  A 


d;spia<cment  n  i  0 )  and  preload.  However,  their  values  are  smaller  than 
those  for  w  =  0.2  for  the  same  value  of  the  freeplay. 
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3.5.2  Effect  Freeplay  for  Constant  Preload 

Figures  27,  22  and  28  show  the  effect  of  varying  the  freeplay 
on  the  limit-cycle  flutter  boundaries  for  a  fixed  preload  of  0.5°  and 
frequency  ratio  u  =  0.8.  Similar  to  the  results  for  u  =  0.2, 
increasing  the  freeplay  will  increase  the  region  where  limited 
amplitude  oscillation  occurs.  It  is  observed  that  the  bulges  on  the 
flutter  boundary  curves  for  a(0)  >  0  diminish  as  the  freeplay  is 
increased.  Also,  the  pockets  of  limited  amplitude  oscillations  which 
occur  inside  the  damped  oscillation  region,  as  found  in  case  5,  are  not 
present  when  w  =  0.8. 

The  limited  amplitudes  a  and  £  are  given  in  Figures  29,  25 

A  A 

and  30  for  values  of  freeplay  0.25°,  0.5°  and  1.0°  respectively. 

Again,  they  are  independent  of  initial  displacement  a(0).  For  a  given 
speed  ratio,  higher  amplitudes  are  obtained  for  larger  values  of 
freeplay.  The  breaks  in  the  curves,  similar  to  those  shown  in  Figures 
13  for  w  =  0.2,  are  not  found  in  the  corresponding  case  for  w  =  0.8 
given  in  Figure  30  because  pockets  of  limited  amplitude  oscillations 
are  not  encountered.  Comparisons  with  results  for  u>  =  0.2  given  in 
Figures  12,  8  and  13  show  that  the  amplitudes  of  oscillation  at 
limit-cycle  flutter  attain  a  larger  value  for  the  smaller  ratio  of  the 
natural  frequencies  u . 

In  section  3.2,  it  is  observed  that  at  u  -  0.2  and  for  a 

particular  combination  of  preload  and  freeplay,  the  results  can  be  used 

for  other  values  of  preload  and  freeplay  as  long  as  the  ratio  of 

preload  to  freeplay  is  the  same.  The  vertical  scale  of  the  flutter 

boundaries  and  limited  amplitudes  a,  and  £ ,  have  to  be  multiplied  bv  a 

A  A 

factor  equal  to  the  ratio  of  the  preloads.  The  results  for  w  =  OH 
show  the  same  behaviour. 
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3.5.  Effect  of  Airfoil-Air  Mass  Ratio  \i 

'Hie  etf»*<  t  of  p  on  the  flutter  boundaries  can  be  obained  from 


and  In  of  Table  1,  and  the  results  are  shown  in  Figures 
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The  limited  amplitude-  n  given  in  Figures  25a,  33a  and  32a 

a  is  slightly  larger  when  u  is  decreased  from  250  to  50. 
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preload  of  0.5°  and  the  values  of  freeplay  being  0.5°  and  1' 


respectively.  The  two  curves  are  identical  for  0.885  <  U*/L’*  <  1.0, 


and  the  curve  for  the  1°  freeplay  extends  to  =  0.835  which  is  the 


location  of  the  limit-cycle  flutter  boundary.  Examination  of  the 
osc  illations  for  these  two  cases  shows  that  the  second  harmonic  is  too 
small  to  be  detected  even  near  the  limit-cycle  flutter  boundary. 

At  zero  preload  with  equilibrium  at  a  point  on  the  linear 
portion  of  tin-  moment -d l sp 1 arement  curve  just  ahead  of  the  freeplay 
case  o  of  Table  1),  Figure  37  shows  the  period  of  oscillations  for 
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For  nonzero  values  of  the  preload,  three  types  of  oscilla¬ 
tions  are  possible,  namely;  damped,  limited  amplitude  and  divergent. 
Divergent  flutter  occurs  when  the  nondimensional  speed  of  the  airflow 


past  the  airfoil  U*  is  greater  than  the  linear  flutter  speed  U*  In 

Li 

spite  of  the  presence  of  a  structural  nonlinearity  of  the  type 
considered,  the  nonlinear  divergent  flutter  speed  is,  within  the 
numerical  accuracy,  equal  to  U*  for  the  cases  studied.  For  U*  less 
than  l-",  limited  amplitude  oscillations  occur  until  a  value  of  U*  is 
reached  where  further  decrease  will  result  in  damped  oscillations.  The 
location  of  the  limit-cycle  flutter  boundary  depends  on  the  airfoil 
parameters  in  addition  to  the  values  of  the  spring  preload  and 
f  r  e  e  p  1  a  y  . 

For  some  combinations  of  preload  and  freeplay,  there  are 
pot.  kets  of  limit-cycle  flutter  in  the  damped  oscillation  regions.  The 
oscillations  inside  these  pockets  are  characterized  by  the  presence  of 
large  sec  on  i  harmonics.  Bringing  the  uncoupled  natural  frequencies  of 
the  two  modes  closer  to  each  other  seem  to  eliminate  these  pockets,  but 
large  bulges  are  detected  on  the  limit-cycle  flutter  boundary.  The 
osc  illations  it.  the  limited  amplitude  region  are  mainly  sinusoidal 
except  vi  ry  near  the  limit-cyc  le  flutter  boundary  where  se  ond 
l.-iT  lifiti  s  are  observed  . 

Th*  limited  amplitudes  and  a  of  the  plunge  and  pitch 

rt  A 

m/t  .on-.  :i  •  .  t  ;  ve  1  v  for  given  ait  foil  parameters  are  dependent  or.  the 
v  .  .1  ■  •  ‘  i  freeplay.  Large  r  am;  1  i  t  tides  are  obtained  win  i  tin  values 
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of  the  freeplay  are  increased.  Changes  in  preload  and  initial 

displacement  do  not  appear  to  have  any  effect  on  them.  The  amplitudes 

also  vary  with  the  ratio  of  the  uncoupled  natural  frequencies  of  the 

plunge  to  pitch  motions  w.  For  the  lower  value  of  w  =  0.2  considered 

in  this  study,  larger  oscillatory  motions  in  both  degrees  of  freedom 

are  observed  than  those  at  w  =  0.8.  Varying  the  airfoil-air  mass  ratio 

p  at  w  =  0.2  does  not  have  any  noticeable  effect  on  a  ,  but  £  is  found 

A  A 

to  increase  with  p.  At  w  =  0.8,  a  decreases  slightly  with  increasing 

A 

p  while  the  corresponding  decrease  in  £  is  much  larger. 

A 

The  period  of  oscillations  at  the  lower  value  of  u  ;  0.2 
decreases  as  U*  is  increased  from  the  limit-cycle  flutter  boundary 
until  a  minimum  is  reached.  Thereafter,  it  increases  at  a  gradual  rate 
towards  the  divergi.it  flutter  boundary.  At  the  higher  value  of 
w  =  0.8,  only  a  gradual  increase  of  the  period  from  the  limit-cycle 
flutter  boundary  to  the  divergent  flutter  boundary  is  observed. 

For  zero  preload,  only  divergent  and  limit-cycle  flutter  are 
observed.  The  divergent  flutter  boundary  occurs  very  close  to  that  for 
the  linear  flutter  case.  Results  for  u  =  0.2  and  p  =  100  show  that 
they  are  independent  of  the  equilibrium  position  of  the  freeplay.  The 
limited  amplitude  curves  versus  speed  ratio  have  a  discontinuity.  For 
speed  ratio  less  than  its  value  at  the  discontinuity,  the  oscillations 
are  characterized  by  the  presence  of  second  harmonics,  while  in  the 
region  between  the  discontinuity  and  the  divergent  flutter  boundary, 
the  osc i nations  are  practically  sinusoidal.  The  variations  of  the 
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period  of  oscillations  with  speed  ratio  in  these  two  regions  are  quite 
different.  A  rapid  increase  of  the  period  with  increasing  speed  ratio 
is  seen  up  to  the  discontinuity.  Thereafter,  the  period  changes  very 
little  as  the  divergent  flutter  boundary  is  approached. 

The  results  presented  in  this  study  on  the  response 
characteristics  of  a  two  degree  of  freedom  nonlinear  system  with 
preload  and  freeplay  are  useful  to  the  understanding  of  the  impact  of 
limit-cycle  oscillations  on  the  structural  integrity  of  aerodynamic 
surfaces.  They  can  also  be  used  to  assess  the  validity  and  accuracies 
of  the  method  of  harmonic  balance  commonly  used  in  the  study  of 
nonlinear  flutter. 
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FIG.  1:  TWO-DEGREE-OF-FREEDOM  AIRFOIL  MOTION 


-f - -4 -  - 4+---  -+ —  — I - 

0  .6  C."7  0.8  '■ .  0 


SPEED  RATIO  U*/U‘ 


FIG.  3:  FLUTTER  BOUNDARY  FOR  m  =  100,  T  =  0.2,  h(<  =  0.25  , 
hf-  =  0.75  AND  Vr  =  0.25 
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FIG  6c  TIME  VARIATION  OF  -  FOR  ..(0)  =  8  ,  U*/Ul  =  0.81 
(CORRESPONDING  TO  LOCATION  3  IN  FIG.  3) 
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FIG  6d  TIME  VARIATION  OF  ■  FOR  .(U)  8  ,U*  Ul  0.82 

(CORRESPONDING  TO  LOCATION  4  IN  FIG  3) 
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FIG.  6e.  TIME  VARIATION  OF  a  FOR  a(0)  =  8  ,  U*/UL  =  0.83 
(CORRESPONDING  TO  LOCATION  5  IN  FIG.  3) 
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FIG.  6f:  TIME  VARIATION  OF  <>  FOR  o(0)  8  ,  U*/ll[  0  95 

(CORRESPONDING  TO  LOCATION  6  IN  FIG  3) 
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FIG.  6g:  TIME  VARIATION  OF  a  FOR  q(0)  =  8  ,  U*/U[  =  0.99 
(CORRESPONDING  TO  LOCATION  7  IN  FIG.  3) 


FIG  6h:  TIME  VARIATION  OF  „  FOR  <.(0)  -  0.3  ,  U*/U‘L  0  99785 

(CORRESPONDING  TO  LOCATION  8  IN  FIG  3) 
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FIG.  7a:  VARIATION  OF  LIMIT  AMPLITUDE  aA  WITH  FIG.  7b:  VARIATION  OF  LIMIT  AMPLITUDE  £A  WITH 

SPEED  RATIO  FOR  m  =  100,  £  =  0.2,  hf+  =  0.25°,  SPEED  RATIO  FOR  m  =  100,  oj  =  0.2,  hf+  =  0.25°, 

hf-  =  0.75°  AND  V,  =  0.25°  hf-  =  0.75  AND  Vf  =  0.25 
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FIG.  11:  FLUTTER  BOUNDARY  FOR  n  =  100,  w  =  0.2,  hf  +  =  0.! 
hf-  =  1.5  AND  Vf  =  0.5 
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12a:  VARIATION  OF  LIMIT  AMPLITUDE  aA  WITH  FIG.  12b:  VARIATION  OF  LIMIT-AMPLITUDE  £A  WITH 

SPEED  RATIO  FOR  n  =  100,  c3  =  0.2,  hf+  =  0.5°  SPEED  RATIO  FOR  v  =  100,  w  =  0.2,  hf+  =  0.5°, 

hf-  =  0.75  AND  V,  =  0.5°  h,-  =  0.75°  AND  V,„  =  0.5° 


13a:  VARIATION  OF  LIMIT  AMPLITUDE  «A  WITH  FIG.  13b:  VARIATION  OF  LIMIT-AMPLITUDE  £A  WITH 

SPEED  RATIO  FOR  n  =  100,  uj  =  0.2,  hf+  =  0.5°,  SPEED  RATIO  FOR  n  =  100,  cj  =  0.2,  hf+  =  0.5°, 


15a:  VARIATION  OF  LIMIT  AMPLITUDE  aA  WITH  FIG.  15b:  VARIATION  OF  LIMIT-AMPLITUDE  £A  WITH 

SPEED  RATIO  FOR  u  =  100,  O  =  0.2,  h,+  =  0°,  SPEED  RATIO  FOR  m  =  100,  C3  =  0.2,  hf+  =  0°, 
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FIG.  18:  FLUTTER  BOUNDARY  FOR  ju  =  250,  -  0.2.  hf  *  0.5  , 

h,_  =  1.0  AND  Vf  0.5 
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FIG  22:  FLUTTER  BOUNDARY  FOR  a  100.  O  B,  h(.  0.5 

h,  1.0  AND  V(  0.5 
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FIG.  23:  FLUTTER  BOUNDARY  FOR  n  =  100,  ~  0.8.  hf .  1.0 

h.  1.5  AND  V,.  1.0 
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FIG.  24a:  VARIATION  OF  LIMIT  AMPLITUDE  aA  WITH  FIG.  24b:  VARIATION  OF  LIMIT-AMPLITUDE  £A  WITH 

SPEED  RATIO  FOR  /j  =  100,  o3  =  0.8,  hf  +  =  0.25\  SPEED  RATIO  FOR  a*  =  100,  a;  =  0.8,  hf+  =  0.25°, 

h,_  =  0.75  AND  V«,  =  0.25  hf-  =  0.75°  AND  Vf„  =  0.25° 
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FIG.  27:  FLUTTER  BOUNDARY  FOR  v  100,  “  0.8,  h, .  0.5, 

hf  0.75  AND  V,  0.5 
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FIG.  28:  FLUTTER  BOUNDARY  FOR  m  100,  0.8,  h, . 

h,  1.5  AND  Vf  0.5 
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ARIATION  OF  LIMIT  AMPLITUDE  <>A  WITH  FIG.  30b:  VARIATION  OF  LIMIT-AMPLITUDE 

RATIO  FOR  n  =  100,  =  0.8,  hf  +  -  0.5  ,  SPEED  RATIO  FOR  n  =  100,  o3  =  0.8,  hf+  = 

h,  1.5  AND  V,  =0  5  h,  =  1.5  ANDV,  =0.5 
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FIG  31  FLUTTER  BOUNDARY  FOR  „  50 
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FIG.  32:  FLUTTER  BOUNDARY  FOR  ij  250,,.:  0  8,  h,.  0  5 
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FIG.  33a:  VARIATION  OF  LIMIT  AMPLITUDE  oA  WITH  FIG.  33b:  VARIATION  OF  LIMIT-AMPLITUDE  £A  WITH 

SPEED  RATIO  FOR  v  =  50.  d  =  0.8,  hf  +  =  0.5  ,  SPEED  RATIO  FOR  n  =  50,  <Z>  =  0.8,  hf+  =  0.5\ 

hf-  =  1.0  AND  V,  =0.5  hf- =  1.0  AND  Vr  =  0.5 


34a  VARIATION  OF  LIMIT  AMPLITUDE  oA  WITH  FIG.  34b:  VARIATION  OF  LIMIT  AMPLITUDE  £A  WITH 

SPEED  RATIO  FOR  M  =  250,  =  0.8,  hf  +  =  0.5  ,  SPEED  RATIO  FOR  n  =  250,  63  =  0.8,  hf+  =  0.5 , 

h,.  -  1.0  AND  V,  -  0.5  h,.  =  1.0'J  AND  V,.  =  0.5 
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The  coefficients  P  ,  P^.  ^21’  ^22  an<^  t^6  terms  and  ^2 


in  Eqs .  (14)  and  (15)  are  given  as  follows: 
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**  ‘  *  V  'j-  ■*-*  O  V  " 


[  #*1  ■**  »,l  wr 


76 


V. 


b: 


is 


yOyfi 


r  ~(i+a.)  (ae  "“‘+ce  idAt)  [  $ "  (t -2Ax  )  +  ( *-a  )a"  (  x -2Ax  )  ♦  a'u-24i)] 
1 2p  ah  n 


-b( x +  Ax )  -d(x+Ax) 


"  r/(K)  (ae  u'l'“w+ce  [V  (0)+(i-ah)a’  (0)+a(0)  ] 


M  a  w  h 


i  ci*+v -~bAT  j  r«i(K> ce 


2  -2, 


-dAx 


(A-9) 


APPENDIX  B 


In  Eq .  (26),  the  elements  of  [Q]  and  [Y  n  }  are: 


p  ( *JXa_ah) 


( B  —  1 ) 


-  ( i+y) 


l  fl  A  2,  •, 

M  ll  8+ah  +  ^ 

r 

a 


(  B— 2  ) 


( B  —  3  ) 


1  -2  ,  , 
-  r  (yx  -a,  ) 
p  a  a  h 


£  [l  +  (i-ah)  [ 2 £  —  2 ( a+c ) ] ]  a'(O) 


[ 2yC^uU*  +  2£-2 ( a+c ) ]  C ' (0) 


(B-4) 


i  [ 2£-(a+c) ]a(0) 


(B-5) 


-  -  i  [2wCaU*_1  +  r"2(|-ah)  [ l-(i+ah) [2£-2(a+c) ] ] ]  a'(0) 

+  J  [r”2^+ah)  [  2£-2(a+c)  ]  ]  C’(0)  +  J  ^2(i+ah)  [2£-2(a+c) 


-  F  (  a  (  0 )  ] 


(B-6) 


[2nCai:*  +  ra  (i~ah)  [l-(i+ah)[2i-2(a+c)]]J  a" 
[r'2(i+ah)  [ 2£-2(a+c)>]  V"  (0) 

[ra2(|+ah)  [ 2£-2(a+c)  +  (|-ah) (2ab+2cd) ] j  a"(0) 

ra2(i+ah)  (2ab+2cd)  5"(0) 

[r"2(|+ah)  [ 2ab+2cd  -  (i~a  ) (2ab2+2cd2) ] ]  o'(0) 

r‘2  (i+ah)(2ab2+2cd2)[C' (0)+a(0)] 
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